FIAN-30/99 
UCRHEP-E255 
14 May 1999 
revised 24 August 1999 



Energy dependence of mean 
multiplicities in gluon and quark jets 
at the next-to-next-to-next-to-leading 

order 

I.M. Dremin 1 and J.W. Gary 2 

1 Lebedev Physical Institute, Moscow 117924, Russia 
2 Department of Physics, University of California, Riverside CA 92521, 

USA 

Abstract 

Analytic predictions for the energy dependence of the mean multi- 
plicities in gluon and quark jets are presented at the next-to-next-to- 
next-to-leading order (3NLO) of perturbative QCD and are compared 
to experiment. The 3NLO correction to the gluon jet multiplicity is 
found to be small. The corresponding theoretical expression provides 
a good description of available gluon jet measurements. The 3NLO 
correction to the quark jet multiplicity is large. The theoretical ex- 
pression for quark jets also describes the data accurately, however, 
but not with the same parameter values as are used for gluon jets. 
It is shown that the well known success of the next-to-leading order 
(NLO) approximation in describing the energy evolution of quark jet 
multiplicity can be attributed to the equivalence of the quark and 
gluon expressions at NLO to within a constant factor, and to almost 
constant contributions from higher order terms to the gluon jet result. 
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Particle multiplicity is one of the most basic properties of a jet. The 
multiplicity of a jet has been the subject of many theoretical and experi- 
mental studies Analytic predictions for the energy dependence of mean 
multiplicities in gluon and quark jets can be derived from perturbative QCD 
(pQCD). The total multiplicity of a jet, i.e. the multiplicity in full phase 
space, is subject to important corrections in pQCD and thus provides a sen- 
sitive means to test higher order perturbation theoryQ In the present study, 
we present a test of higher order perturbative calculations based on the total 
multiplicity in jets. 

The perturbative solutions for multiplicity are usually expressed in terms 
of the anomalous dimension 7 and the multiplicity ratio r, defined by 

i(y) = — rr= ln ™G(y) ; ny) = — rr > i 1 ) 

now n F {y) 

where hq and n-p are the mean total multiplicities in gluon and quark jets, 
respectively. In these expressions, y=ln (jo#/Q ) defines the energy scale of 
the jets, with p the momentum of the parton which initiates the jet, 9 the 
opening angle of the first branching in the jet, and Qo a cutoff which defines 
the limit of perturbative evolution. For p large and 9 small, the virtuality 
Qjet (invariant mass) of the jet is given by Qj e t~P#/2. Primed quantities 
in (0) denote derivatives with respect to y: n' G =driG/dy, etc. 

Perturbative expressions for uq and n-p at the next-to-leading order (NLO) 
and the next-to-next-to-leading order (NNLO) approximations of pQCD have 
been presented in || and respectively. Here we present the solutions at 
the next-to-next-to-next-to-leading order (3NLO) approximation and com- 
pare the results to experiment. 

The perturbative expressions for 7 and r at 3NLO are: 

j(y) = 7o(l - ai7o - a 2 7o - a37o) + C(7o) (2) 
r{y) = r (l-r 1 7 -r 2 7 2 -r 3 7o 3 )+0(7o 4 ) , (3) 
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-'^This is in contrast to multiplicity in a limited region of phase space, such as the 
multiplicity of soft particles, which is less sensitive to higher order corrections; this fact 
has recently been exploited to determine the color factor ratio Ca/Cf using the soft 
particle multiplicities in gluon and quark jets || (see also 
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where as is the strong coupling strength, /3o=(llNc — 2Np)/3, /3i=(51Nc — 
19Np)/3, Nc=3 is the number of colors, Np is the number of active quark 
flavors, r =N c /CF, and Cp=4/3. The coefficients and are calculable in 
pQCD and are given in || for 2=1,2,3.0 

The expressions for 7(2/) in ([I]) and (^) can be combined to yield 



n G (y) oc exp 



y 
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(5) 



In conjunction with (Q), we then determine the energy dependence of the 
gluon jet multiplicity at 3NL0 to be: 



n G {y)=Ky- aiC 'exp[2C > /y + S G {y)] , (6) 



where K is an overall normalization constant, C=-i/4Nc//9o ; an d 



^(2a 2 C 2 + §[ln(2y) + 2] 




1 M2y) + 1] • (7) 



To obtain this result, the term proportional to f3\ in parenthesis in equa- 
tion (£|) is treated as being small compared to 1. The first term in the 
exponent in equation (|6|) is the leading order (LO) term. The term in front 
of the exponential arises at NLO. The term proportional to in flj]) 

enters at NNLO, while that proportional to 1/y is the 3NLO contribution. 

Experimental measurements of the inclusive charged particle multiplicity 
of gluon jets as a function of energy scale Q are shown in Fig. |l| These results 
utilize a definition of gluon jets which corresponds to that employed for the 
calculation, based on the production of a virtual gluon jet pair, gg, from a 
color singlet point source. These are the only such gluon jet data currently 
available. The three data points at scale Q ~ 5 GeV are derived from the 
hadronic component of T(1S)— > 7gg decays J7J. The virtuality Q is given by 
the invariant mass of the hadronic system. Similarly, T(3S)— > 7Xb2(10.27)^ 
7gg decays provide the measurement at Q ~ 10 GeV [§] , with the scale given 
by the xvi mass. The result at Q ~ 80 GeV p[ is based on hadronic Z° data: 

2 For Np=3, the value we emphasize in this letter, 01, a-i and 0,3 are 0.280, -0.379 and 
0.209, while n, r 2 and r 3 are 0.185, 0.426 and 0.189 fj|. 
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Z° — > qqginci., in which gi nc i. refers to a gluon jet hemisphere recoiling against 
two almost collinear quark jets q and q in the opposite hemisphere. The g inc i. 
jets correspond closely to hemispheres of color singlet gg events, as discussed 
in prof] . For the results shown here, the gi nc i. hemisphere results reported 
in have been multiplied by a factor of two both for the multiplicity and 
energy scale so that they correspond to full gg events analogous to the T 
data. 

The solid curve in Fig. [l] shows the result of a \ 2 fit °f expression (^) 
to the data, with j/=ln(Q/A) where A is a fitted parameter. The other fit- 
ted parameter is the normalization constant K in (||). For this fit, Nf=3. 
The T(1S) measurements J7| are not included in this fit because a system- 
atic uncertainty was not provided for them. Nonetheless, these points are 
seen to lie near the fitted curve. The results for the fitted parameters are 
A=1.03 ± 0.24 GeV and K=0.288 ± 0.037. The uncertainties are defined by 
the maximum deviations observed when the gluon jet measurements are var- 
ied by their one standard deviation uncertainties. These results, along with 
those found using Np=4 and 5, are summarized in the top portion of Table |l[ 
The results of the fits with Np=4 and 5 are virtually indistinguishable from 
those shown by the solid curve in Fig. [1|. Note that the result using Np=3 
probably has the most physical relevance since b and c quarks are rarely 
produced in the perturbative evolution of the jets, even for the Z° data. The 
dashed curve in Fig. |T] shows the prediction of the Herwig parton shower mul- 
tihadronic event generator [ITf, version 5.9, for the inclusive charged particle 



multiplicity of gg events as a function of the center-of-mass energy, Q=E C 
The Monte Carlo parameters are the same as those used in 0. Herwig is 
seen to describe the data well and to be similar to the 3NLO result. 

The analytic expression for quark jet multiplicity can be obtained from (|l|) 
and©: 

/ \ n G(y) n G (y) 

"Flf = = — n 2 sT ' ( 8 ) 

with noiy) given by (|6]). The expression for r(y) in the denominator of (§) 
is known to give a relatively poor description of data, however. For example, 
the 3NLO prediction at the scale of the Z° is r ~ 1.7 ||, about 13% larger than 
the experimental result of 1.51 ±0.04 0. Therefore, it can be anticipated that 
the description of quark jet multiplicity provided by (§) will be deficient if 
the values of A and K found using the gluon jet data are employed. It should 
be noted that this problem with r is even worse at lower orders, e.g. the 
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NLO and NNLO predictions for r are about 2.1 and 1.8, respectively, in 
even greater disagreement with experiment than the 3NLO result. 

Experimental measurements of quark jet multiplicity are shown in Fig. |2] 
[|T2]]. These data are the inclusive charged particle multiplicity values of 
e + e~— > (Z°/7)* — > hadrons events and correspond to the definition of quark 
jets employed for the calculation, i.e. the production of a virtual qq pair 
from a color singlet. The scale is Q=E cm .. The results shown for the LEP 
experiments are combined values of ALEPH, DELPHI, L3 and OPAL. The 
combined values are obtained using the unconstrained averaging procedure 
described in [[n| , for which a common systematic uncertainty is defined by the 
unweighted mean of the systematic uncertainties quoted by the experiments. 
LEP-1 refers to data collected at the Z° peak, LEP-1.5 to data collected at 
E c . m .~133 GeV, and LEP-2 to data collected at or above the threshold for 
W + W~ production. 

The vertically striped band in Fig. |2] shows the prediction of (§) for Np=3 
using the values of A and K from the fit to the gluon jet measurements 
(Fig. [I]). The width of the band corresponds to the uncertainties in the 
values of A and K presented above. Almost identical results to those shown 
by the band are obtained using N F =4 or 5. The band typically lies 15-20% 
below the data mostly as a consequence of the problem with the theoretical 
prediction for r noted above. 

It is also of interest to fit expression (|8]) to the quark jet measurements. 
The results of such a fit using Np=3 are shown by the solid curve in Fig. ^. 
The values obtained for A and K are 0.35 ± 0.12 GeV and 0.222 ± 0.030, 
respectively, where the uncertainties are defined by the maximum difference 
between the results of the standard fit and those found by fitting only data 
between 29 and 189 GeV, between 12 and 161 GeV, or by excluding the 
LEP-1 data point. The fit is seen to yield a good description of the data. 
The results found using Np=4 or 5 are essentially identical to those shown by 
the solid curve in Fig. || The values of the fitted parameters are summarized 
in the central portion of Table [1]. The values of A obtained from fitting the 
quark jet data are seen to be three to four times smaller than those obtained 
from fitting the gluon jet data. 

The dashed curve in Fig. |2| shows the prediction of Herwig for the inclusive 
charged particle multiplicity of hadronic events in e + e~ annihilations. The 
Herwig prediction reproduces the measurements well, analogously to what 
was observed for gluon jets. 

The expression for n F given by equation (H) is not entirely satisfactory 
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from a theoretical perspective. The purely perturbative 3NLO result for rip, 
analogous to equation (|G|) for uq, is 



n F (y) = ^-y- a ^ 2 exp[2C^ + 6 F (y)} , (9) 



with 

S F (y) 



|r 1 + 2a 2 C 2 + ^[ln(2y) + 2]| 



C 

+y{a3C 2 + ^ + r 2 -^[ln(2y) + l]} . (10) 

where the same notation is used as in equations @ and (|7]). Equation (|]) 
is derived from equation (|), keeping terms only up to y~ l in the exponent. 
It is seen that a term proportional to r 3 appears in ([Bp but not in (|9|): such 
a term would be proportional to y~ 3 ^ 2 in (10) and thus would contribute to 



the 4NLO approximation. This difference between (|8|) and (|9|) arises because 
of the difference in the order of 70 in (0) and (0). The manner in which the 
powers of 70 are structured in (0) and ([|) is conventional because 7 enters in 
the exponent whereas r is used as a multiplicative factor. As a consequence 
of this difference between (|2]) and (|3]), all corrections of a given order in 
the denominator of (^) contribute to yet higher orders in ([|). For example, 
r% contributes only to the NLO term in the denominator of (|8[) whereas it 
contributes to the NNLO (~ y _1 / 2 )and 3NLO (~ y~ l ) terms in (§). This 
situation can also be considered to arise from the fact that the logarithmic 
slope of quark jets, j-p, is given by 



7f 



[lnn F ]' = 7 -- • (11) 



r 



However, from (|3|) and (|J), r' ~ — r^ri^ ~ To r i7o- Therefore, the first 
correction to 7f appears at NNLO, but it is determined by the NLO correc- 
tion to r. We conclude that it is improper to use the term with r 3 in @ 
until the 4NLO contribution to 7 is known. By extension this implies that 
if the NLO formula is used to describe quark jet multiplicity, as is common 
practice (see below), the LO result r=ro should be inserted into all formulas 
to be self-consistent within the perturbative approach. Therefore the mean 
multiplicities in quark and gluon jets coincide in NLO up to the constant 
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normalization factor ro=9/4. Any usage of r different from r$ implies that 
NNLO corrections have been included. 

The dotted curve in Fig. |2] shows the result of a x 2 fit of equation (Q) to 
the quark jet data for N F =3. The latter equation is seen to provide a good 
description of the measurements, equivalent to that given by equation 
The fitted curves for Np=4 and 5 are indistinguishable from the curve for 
Np=3. The values of the fitted parameters are given in the bottom section 
of Table |l[ These values are similar to those found from fitting equation (||) 
to the quark jet data. 

Our method of comparing the 3NLO predictions to data is in contrast to 
the procedure usually employed to test QCD analytic predictions for multi- 
plicity. Typically, the NLO expression for gluon jets |4| is used to fit quark 
jet data. A fit to gluon jet data is rarely performed. At NLO, the theo- 
retical expressions for rip and hq differ by only the constant factor r , as 
noted above, which justifies this procedure. The NLO expression is found to 
provide a good description of the measurements, analogous to the agreement 
we find in the fit of equations (||) or (|9|) to quark jets (the solid or dotted 
curves in Fig. ^j. This good description belies the fact that the appropriate 
prediction for r at this order, ro=2.25, is in striking disagreement with the 
experimental value r^l.5. 

The good description of quark jet multiplicity at NLO is easily under- 
stood, because the sum of the NNLO and 3NLO corrections to n G , given 
by Sg{v) (equation (^)), is small and almost constant at presently accessible 
energies. To illustrate this, we show in Fig. |3] the behavior of £g(z/) as a 
function of scale for Np=3, 4 and 5. For the results of this figure, the scale is 
2/=ln(Q/A) for which the values of A are taken from the fits to quark jet data 
using equation ([|]) (central section of Table [l]) to correspond to the procedure 
outlined above to test the NLO equation. It is seen that Sq is almost constant 
for the range of energies relevant to experiments. Furthermore, 5q is small 
compared to the LO term in the exponent of (||), which typically has a value 
from 5 to 7 over the range of scale Q shown. Therefore the NLO formula for 
hq is not modified significantly by the NNLO and 3NLO contributions and 
the expression n F =n G /r fits the quark jet data well irrespective of whether 
the higher order terms are included in uq. 

In summary, we have presented analytic expressions for gluon and quark 
jet multiplicity as a function of energy scale at the next-to-next-to-next-to- 
leading order (3NLO) of perturbation theory. We have performed fits of these 
expressions to gluon and quark jet measurements, using data for which the 
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experimental and theoretical definitions of jets coincide. To our knowledge, 
this is the first time such a fit has been performed for gluon jets. We find 
that the gluon jet expression describes the available gluon jet measurements 
accurately, with values of the parameter A in the range from about 0.6 to 
1.0 GeV depending on the number of active quark flavors, Np. The 3NLO 
expression for quark jets also provides an accurate description of the data 
with values of A about three times smaller than those found for gluon jets. 
We attribute the well known agreement of the next-to-leading order (NLO) 
expression for quark jet multiplicity with data to the facts that the gluon 
and quark jet predictions are equivalent at this order to within the constant 
factor r =9/4, and that the gluon jet prediction is only slightly modified by 
the NNLO and 3NLO terms. 
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A (GeV) K xVd.o.f. 


(a) Gluon jets, equation (EJ) 


N F =3 
N F =4 
N F =5 


1.03 ±0.24 0.288 ±0.037 0.01/2 
0.84 ±0.21 0.244 ±0.034 0.01/2 
0.64 ±0.17 0.205 ±0.031 0.01/2 


(b) Quark jets, equation fl8|) 


N F =3 
N F =4 
N F =5 


0.35 ±0.12 0.222 ±0.030 5.8/13 
0.233 ±0.086 0.170 ±0.024 5.9/13 
0.135 ±0.055 0.126 ±0.019 6.1/13 


(b) Quark jets, equation (^) 


N F =3 
N F =4 
N F =5 


0.322 ±0.098 0.216 ±0.025 5.6/13 
0.231 ±0.083 0.169 ±0.024 5.8/13 
0.148 ±0.059 0.129 ±0.020 6.1/13 



Table 1: Results of a two parameter fit of QCD expressions for the scale 
evolution of event multiplicity to the measured mean charged particle multi- 
plicities of (a) gg (top), and (b) and (c) qq events (center and bottom) from 
a color singlet point source. For the gluon jets, the QCD expression used is 
the exact 3NLO result for n G , equation (Bp. For the quark jets, the QCD 
expression is given either by (b) n-p=nQ/r, equation (§), which contains some 
terms beyond 3NLO as discussed in the text, or by (c) the exact 3NLO result 
for n F , equation (|). A is an energy scale parameter related to the pertur- 
bative cutoff Qo while K is an overall normalization constant. For the gluon 
jets, the uncertainties of A and K are evaluated by varying the gluon jet data 
points by ± 1 standard deviation of their total uncertainties. For the quark 
jets, the uncertainties are evaluated by varying the fit range as described in 
the text. 
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20 40 60 80 100 120 

Scale Q (GeV) 

Figure 1: The mean charged particle multiplicity of gg events from a color 
singlet point source versus energy scale Q. The solid curve shows a fit of the 
3NLO expression to the data using Np=3, where Np is the number of active 
quark flavors. The dashed curve shows the prediction of the Herwig Monte 
Carlo for gg events. 
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Figure 2: The mean charged particle multiplicity of e + e~ hadronic annihi- 
lation events versus energy scale Q=E c m . The dotted curve shows a fit 
of the 3NLO expression to the data using Np=3. The hatched band and 
solid curve show results using the expression nF=nc/r as explained in the 
text. The dashed curve shows the prediction of the Herwig Monte Carlo for 
e + e~ — > qq events. 
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Figure 3: The behavior of the sum of the NNLO and 3NLO correction term 
to the gluon jet multiplicity, 5q, as a function of energy scale Q, for Nf=3, 
4 and 5. 
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